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We study the bosonic atoms with a wide Feshbach resonance at zero temperature in terms of the 
renormahzation group. We indicate that this system will always collapse in the dilute limit. On the 
side with a positive scattering length, the atomic superfiuid is an unstable local minimum in the 
dilute limit and it determines the thermodynamics of this system within its lifetime. We calculate 
the equilibrium properties at zero temperature in the unitary regime. They exhibit universal scaling 
forms in the dilute limit due to the presence of a nontrivial zero temperature, zero density fixed 
point. Moreover, we find that the T — thermodynamics of this system in the unitary limit is 
exactly identical to the one for an ideal Fermi gas. 
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I. INTRODUCTION 

The bosonic atoms with Feshbach resonances have 
been one of the exciting ultrocold systems studied so 
far. On the side with a negative scattering length, the 
controlled collapsing and exploding dynamics has been 
observed^. On the other side where the scattering length 
becomes positive, molecules can be produced and this 
system becomes a mixture of atoms and molecules. De- 
spite the reduced lifetime of molecules in the unitary 
regime due to the enhanced three-body recombination^, 
quite remarkable progress has been made and the prop- 
erties of the strongly repulsive ultracold bosonic atoms 
have also been studied experimentally^. 

A mean-field phase diagram for the Bose gas with an 
s-wave Feshbach resonance has been proposed^. At zero 
temperature, there are two phases: the atomic super- 
fluid (ASF) in which both atoms and molecules condense 
and the molecular supcrfluid (MSF) in which only the 
molecules condense. Hence, there is a quantum phase 
transition between the ASF and MSF phases^i^. The 
stability of these phases at T = was examined within 
the mean-field level^. Recently, a Nozieres-Schmitt-Rink 
formalism was developed to study the finite tempera- 
ture properties of this system^. The ground state of 
the cold Bose gas with a large positive scattering length 
has also been studied by the variational wavefunction 
approach^. Near the Feshbach resonance, various phys- 
ical quantities were found to exhibit universal scaling 
forms^. 

Although the idea of the Feshbach resonance has its 
root in the few-body physics, it was shown that, for a 
dilute quantum gas, it is possible to cast the low-energy 
physics near the Feshbach resonance and the associated 
universal properties into the framework of the renormal- 
ization group (RG) due to the presence of a nontrivial 
zero-density fixed pointi^. Bearing this picture in mind, 
the low-energy properties of a dilute Fermi gas with a 
wide Feshbach resonance in the unitary regime can be 
computed either by an e expansionii (e = 4 — c?) or a 
large N expansio n^°d^ . Here we will extend the idea in 



Ref. [lOl to the dilute Bose gas with a wide Feshbach 
resonance. In contrast to the case of fcrmionic atoms, 
the self-interactions between atoms and molecules must 
be included in order to guarantee the stability of the 
uniform phases although these operators are irrelevant 
in the sense of RG. This point was also emphasized in 
Ref. [^. The resulting stability conditions impose a lower 
bound on the value of the chemical potential for atoms. 
Hence, we conclude that a Feshbach-resonant (FR) Bose 
gas will always collapse in the dilute limit. Moreover, in 
the dilute limit, the ASF state can appear in the guise 
of an unstable local minimum for the free energy den- 
sity, while the MSF state cannot exist. The ASF state 
will eventually collapse. However, within its lifetime, the 
properties of this system may be captured by the ASF 
state when we approach the unitary limit from the side 
with a positive scattering length. Then, we use the e ex- 
pansion to calculate various properties of the ASF state 
at T = 0, including the relation between the chemical 
potential and detuning, the equation of state, and the 
sound velocity. We show that the thermodynamical prop- 
erties at T = exhibit universal scaling forms near the 
Feshbach resonance due to the presence of a nontrivial 
zero-temperature, zero-density fixed point. Surprisingly, 
we find that the T = thermodynamics of the ASF state 
in the unitary limit is exactly identical to the one for an 
ideal Fermi gas. This correspondence was also noticed 
m Refs. i and i from a different approach. Our main 
results are summarized in Figs. [Tl[2l[3l andlH 



The rest of this paper is organized as follows: In Sec. 
im we present the RG theory to fix the notation. We 
discuss the stability of the dilute Bose gas with a wide 
Feshbach resonance and its phase diagram at T = in 
Sec. mil The equilibrium properties of the ASF state 
at T = are calculated in Sec. IIVI The last section is 
devoted to a conclusive discussion and comparison with 
previous works. 
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II. THE RENORMALIZATION-GROUP 
THEORY 

We start with the action for the atom-molecule model 
in the grand-canonical ensemble^: 



dr d'^x U 



(1) 



where 



U = .go + H.c.) +^/3|*a|'|*. 



Here '9 a a-nd are the annihilation operators of atoms 
and molecules, respectively, m and M are the masses 
of the atom and the molecule, respectively. Galilean in- 
variance requires that M = 2m. We will see later that 
So = does not correspond to the unitary limit (a j^ — > 
where as is the scattering length between atoms) due to 
renormalization. Hence, we define 



U2 I 



where 5c is a countertcrm which is chosen such that 5 = 
just corresponds to the unitary limit. We further define 

fim = 2)1- S. 

We see that the action S [Eq. ([T])] has a quantum 
critical point (QCP) at /i = = (5 and T = f//3 = 0, 
which will be dubbed as the Gaussian fixed point. That 
is, S* is invariant against the scaling transformation 



xe 



re 



„d//2 



(2) 



at this point when we take go ^ = with a = f , 2, 3. 
We will employ this Gaussian fixed point as a point of 
departure to study the low-energy physics of S in the 
dilute limit and in the unitary regime. That is, we will 
work under the assumption 



H,l^l« 



1 



(3) 



where ao is the background scattering length for atoms. 
We will take ao as the short-distance cutoff for the RG 
flow. 

Now we would like to calculate the RG equations at 
T ~ 0. Under the scaling transformation ([2]), ^, fJ.m, go, 
and Ua transform like 



fi — fie 



Mr, 



21 



9o = 9oe 



el/2 



and 



(e-2); 



where e = A — d. (5c = at the tree level.) Hence, the 
coupling constants UaS are irrelevant as long as d > 2. 
However, to guarantee that the Hamiltonian of this sys- 
tem is bounded from below, we cannot simply set Mq, = 0. 
To stabilize this system, it suffices to keep the ui and U2 
terms. In the following, we will ignore the U3 terrfti^. 

To proceed, we have to calculate the one-loop RG equa- 
tions for /X, fim, and go. Within the e expansion, we find 



dri 
~dl 



2ri 



dr2 „ , ^^''2 



-A 



A3 



2(l-ri)2 



(4) 
(5) 
(6) 



with the initial values 



niO) = 2mfial , r2(0) 
A(0) - y/2IQmgoaf\ 



2m ^, 



m "-0 ' 



where Kd = 2/[{ATTY/'^T{d/2)]. At this order, 5^ is 

d'^p 



Scaling stops when I = where is defined by |7'i(Z*)| 
1. From the solution of Eq. ([4]), we get 



^j2m\fjL\dl 



(7) 



The diluteness condition (the first inequality in Eq. dH)) 
leads to e'* » 1. Since before scaling stops \r\{l)\ <^ 1, 
Eqs. ([5|) and ^ can be approximated as 



-f ^ (2- A2)r2+2A2ri , 



dr2 

~dl 
dX 
H 



-A 



A3 
2 



(8) 
(9) 



We want to emphasize that his approximation does not 
change the fixed-point structure of Eqs. ([4]) - ([6]). In the 
following, we will work with Eqs. dH), ([5]), and (O. 

For d < 4, Eqs. dH), and (O have two fixed 
points: the Gaussian fixed point (rQ,A) = (0,0) and 
the interacting zero-density fixed point {ra,X) = (0,A*) 
with A* = ^/e. The Gaussian fixed point is IR unstable, 
while the interacting zero-density fixed point is IR stable. 
Hence, the latter controls the low-energy physics of the 
dilute FR Bose gas in the unitary regime. 

To relate 6 to the physical detuning ly = — l/a^., wc will 
calculate the T-matrix at zero density, and the result is 



T-\k) 



fc2 - 2m{6 + 6c) 
Amgl 



+Kdm / dq 



q- 



- A:2 _ iO+ 
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Using T(0) = ATias/m and the one-loop result for 5c, we 
get 



47ra, 



mv 



S 

Wo 



By definition, the low-energy physics of the system with a 
wide Feshbach resonance at zero density is parametrized 
by only one parameter a^. Hence, for wide Feshbach 
resonances, we may simply take A(0) = A* if we are only 
interested in the physics in the unitary regime^^, and 
thus we obtain 



(10) 



bmao 



in d = 3 (e = 1) where b = AirKa/Xl = 2/{t:XI). Since 
A(/) = A* for wide Feshbach resonances, Eq. ([8]) can now 
be solved easily. Using Eq. ([7]) to eliminate , we get 

r2iL) = 2sgn(Ai) - (2ma2)'/25|^|i-'^/2 , (11) 

where sgn(a;) = 1,-1 for a; > and a; < 0, respectively. 



III. THE PHASE DIAGRAM AT ZERO 
TEMPERATURE 



The free energy density / is given by 
/ 



-(d+2)i. 

-f{Q , 



2ma 



d+2 



(12) 



where f{l*) is the free energy density for the renormal- 
ized Hamiltonian at / = Since A(Z») = 0(e) and 
Uaih) ~ 0, we may calculate f{h) in terms of the per- 
turbation theory in the coupling constants. The leading 
order results will be given by the mean-field theory of 
the rcnormalizcd Hamiltonian. Moreover, the prefactor 
in front of f{l*) in Eq. ()12|1 is always positive, the phases 
of this system can be determined directly from f{h). 



A. Stability conditions 

To proceed, we make a change of variables: 

X = aox , T = 2malf , ^a,m = a'o'^^^'^ a,m , 
the action S [Eq. ([T])] can be written as 



S 







where /3 = /?/(2ma2), si.2(0) 



and 



At large values of Iv^a,™!, the rcnormalizcd potential U 
a.t I = U is of the form 

16s% 



Kd 



where s* = Sq (/,). In order that the renormalized Hamil- 
tonian is bounded from below, we must require that 

> 0. 

The one-loop RG equations for si^2 are given by 



A^ 



(13) 
(14) 



In Eq. (jl4p . we have neglected the ri and r2 dependence, 
and thus they are valid only in the dilute limit and in 
the unitary regime. Moreover, we do not consider the 
loop corrections due to the ui^2 terms because they are 
irrelevant operators. For the wide Feshbach resonances, 
we may simply set A = A, , and thus we get 

Kd-mui 



1 c d.-2 

loag 

Kdmu2 

loan 



d-2 



(15) 



in the dilute limit. 

The stability conditions require that Mi > and U2 > 
Ur where 



16e2 



(2m|/i|)i-^/2 



{d — 2)Kdm 

In other words, U2 must be strong enough to guarantee 
the global stability of this system in the thermodynamical 
limit. For given U2 > 0, the constraint U2 > Uc can also 
be expressed as a lower bound for i.e. 

2/(d-2) 

(16) 



27Ti|^|ao > 



16e2ag-2 



{d - 2)Kdmu2 



Equation ()16p suggests that the FR Bose gas is stable 
either at moderate density or with strong dimer-dimer 
repulsions. 

The point that the molecule-molecule repulsion is cru- 
cial to stabilize this system was also emphasized by the 
previous mean-field stability analysis^. In the /i, (5 < re- 
gion, the mean-field theory predicts that the ASF state 
is stable for U2 > 0^. Our RG analysis, however, indi- 
cates that in this region this system is thcrmodynam- 
ically unstable as \fi\ smaller than some critical value. 
We notice that this conclusion is consistent with an ear- 
lier RPA calculatioriF^. Within our approach, the reason 
behind this instability is that the molecule-molecule in- 
teraction becomes attractive at low energy even when the 
bare value of U2 is positive. 
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B. The mean-field theory at I — I, 

To calculate f{h), we perform a mean-field theory on 
the action S at I = l^. By inserting the mean-field ansatz 
^a,m('^^j i) — ^a,m into S and then taking the variations 
with respect to $jj , we obtain the mean-field equations 

-rr$, + 2,/jA,$t$„ + ^|$„|2$^ = 0, 

V J^d 



if/ " Kd 



^ml^^ra = ,(17) 



where 2 = fi,2(^*)- We will solve Eq. (fTT)) within 
the e expansion. From Eq. (jlSp . we may neglect the 
si term and the term is necessary only for the ex- 
istence of the MSF state. Hence, there are three types 
of solutions for Eq. (fTT]) : the ASF state characterized 
by |$,| = ^KdrM/i2\,) and |$™| = ^^IQ/8rl/X,, 
the MSF state characterized by $a = and = 
^/Kdrf/(32s^ , and the vacuum state (denoted by N) 
characterized by $a = = The solution for the 

ASF state can exist only for r* 2 > 0, while the solution 
for the MSF state can exist only for 7'2 > 0. In terms of 
these solutions, we find that fN{l*) ~ and 



fASFih) = - 



Kdr'2 



fMSF{l*) 



64sS 



(18) 



The phase diagram can be determined by comparing 
f{U) for the ASF, MSF, and N states. A straightforward 
calculation gives rise to the following results: (i) The ASF 
state occupies the region with rjf > and < < 853 /e; 
(ii) the MSF state occupies the region with > 853 /e > 
0; and (iii) the rest in the r^-rj plane is the N state. 
Substituting the expressions for and [Eqs. (fTT|) 
and (jlSp ] into result (i), (ii), and (iii), and taking into 
account the stability condition [Eq. (fT6|) ]. we obtain the 
phase diagram at T = 0. The ASF phase exists in the 
region with 



2m^aQ > 



{d - 2)Kdmu2aQ 



2-d 



(19) 



and 



[2 - Cd(2TO^ao)](2m/iao) 2 ^ < 2mJaQ < 2{2m^a^)i' 

The MSF phase exists in the region with 

2m5al < [2 - Cd{2mnal)\{2m^jialf/'^^'^ , 



(20) 



for 2m/^aQ > 



and 



2m8al < -2{2m\fi\alf^^~ 



(21) 



(22) 



for 2m/iaQ < — 



Cdix) 



Kdmu2al '' 1 |d/2- 



Here 







1 " 










/ ASF 








MSF \ 


2m/jal 


MSF 











2e 



d-2 



FIG. 1: (color online) The schematic phase diagram of a dilute 
FR Bose gas in d = 3. The unstable region lies between the 
two dashed lines. In the shaded region, the ASF state is an 
unstable local minimum of the free energy density. 



The rest of the phase diagram is occupied by the vacuum 
state. 

A schematic phase diagram is shown in Fig. [1] The 
MSF phase can exist only at finite density of atoms. Fur- 
thermore, its thermodynamical properties do not exhibit 
universal scaling forms even in the unitary limit {6 = 0) 
due to the presence of a dangerously irrelevant operator 
U2- This can be seen from /j\/5f(Z*) [Eq. Although 
the thermodynamical instability in the region with fi > 
and 2mSal < 2{2mfial)'^/^^^ (the shaded region in Fig. 
[1]) implies that the ASF state will eventually decay, the 
time scale, which is governed by the kinetics and dissi- 
pation, may be long enough that the ASF state appears 
stable. The detailed mechanism for the decay of the ASF 
state, however, is beyond the present study. It is hoped 
that the ASF state can survive for a finite lifetime if we 
prepare the system on the side with a positive scattering 
length. (After all, the experiments on ultracold atoms 
involve states which are thermodynamically unstable be- 
cause the ground state of alkali atoms at nano-Kelvin 
temperatures is in fact a solid.) Finally, we must empha- 
size that the exact location of the unstable region (the 
region between the dashed lines in Fig. [1]) is sensitive to 
the short-distance physics, and the determination of it is 
beyond the scope of the present approach. 



IV. THE EQUILIBRIUM PROPERTIES AT 
ZERO TEMPERATURE 



From the above analysis, the dilute FR Bose gas may 
be described by the ASF state on the side with a positive 
scattering length though the ASF state will eventually 
collapse. Hence, we will explore the thermodynamics of 
the ASF at T = 0. 

Inserting fASPih) [Eq. (fT5|) ] into Eq. (fT^ and using 
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FIG. 2: (color online) The chemical potential fj, (in units of FIG. 3: (color online) The pressure P (in units oin^^^/{2m)) 
•n?^^ /{2m)) as a function of ;^ = — 1/os (in units of v}^^) for as a function oi v = — 1/os (in units of for given atom 

given atom density n. density n. 



Eq. (fTTj). we get the free energy density for the ASF state 



Iasf 



- — (2m) 



(23) 



From the thcrmodynamical relation n = —{df/dfi), we 
find that 



1 



2m5a% 



(2to^)^ 



3-e/2 

where n is the density of (bare) atoms and 
Kd{i - e/2) 



(24) 



4e 



Using Euler's relation / = — P, we get the equation of 
state at T = for the ASF phase 



4e 



1 - 



(2m//)'*/2 



(25) 



with ^ as a function of n and 5 given by Eq. ((24|) . Fi- 
nally, we may calculate the sound velocity through the 
relation v1 — (n/m)dfi/dn, and obtain 



2fi 
dm 



{3- e/2){2mny/^-^ -2mSal 
(3 - e/2)(2m^)'^/2-i - imSaf^/d 



(26) 



with ^ as a function of n and (5 given by Eq. (j24p . 

The atom density n gives rise to a characteristic length 
scale n~'^/'^ and a characteristic energy scale r? I j {2m) . 
By measuring the lengths and energies in units of these 
characteristic scales, setting c? = 3, and using Eq. ((TO)) to 
replace 6 by Eqs. p4|) . (|25|) . and (|26|) can be written 
as 



167r2 



,,3/2 



27r 



-xy 



(27) 




FIG. 4: (color online) The sound velocity Vs (in units of 
vi}^'^ l{2ni)) as a function oi v — ~l/aa (in units of n^^^) 
for given atom density n. 



where x = bmSaon = n ^^'^v with v ~ —1/as, b = 
2/tt and y = 2miin~^^^ , and 



p = 


1 

8^ 




TT 




42/ 




2'KX 




- Attx ^ 



(28) 
(29) 



where P = 2mn~^^^P and Vs ~ 2mvsn~^/^ . Equations 
([27|) — p9p exhibit the universal scaling forms for /z, 
P and Wg. These universal scaling forms are due to the 
presence of a nontrivial zero-density fixed point. The 
results are shown in Figs. [2l [3l andlH 

In the unitary limit (5 = Eqs. (f24|l and ((25|) reduce 

to 



1 / n 

2m \Cd 



2/d 



(30) 



6 



„l+2/d 

P = 27d ■ (31) 

We notice that Eqs. ([50]) and (PT|) are identical to the 
Fermi energy and the equation of state at T = for a 
free Fermi gas with mass m* = [dCd/ KdY^'^m. That is, 
within the e expansion, we have shown that the thermo- 
dynamics at r = for a unitary dilute Bosc gas is identi- 
cal to the one for a free Fermi gas. This correspondence 
between the thermodynamics at T = for the dilute FR 
Bose gas and the free Fermi gas was also noticed in Rcfs. 
d andl^ from an entirely different approach. Therefore, 
such a "fermionization" phenomenon should not be an 
artifact of either approach, and deserves further investi- 
gations. 

V. CONCLUSIONS AND DISCUSSIONS 

We study the ground-state properties of a dilute FR 
Bose gas from the point of view that this system can be 
described by a nontrivial zero-density fixed point. From 
a RG analysis, we first show that the FR Bose gas will 
always collapse in the dilute limit. Following from this 
consequence, the MSF state cannot exist in the dilute 
limit while the ASF state can survive with a finite life- 
time if we prepare this system from the side with a pos- 
itive scattering length. Hence, the low-energy properties 
of a dilute FR Base gas in the unitary regime are sup- 
posed to be described by the ASF state. Based on this 
observation, we show that the equilibrium properties of 
a dilute FR Bose gas at T = exhibit universal scaling 
forms near the Feshbach resonance. 

Most of previous works were concentrated on the 
mean-field phase diagram. In the present paper, wc focus 
on the possible universal properties of this system near 
the unitary limit. This point was less studied before, and 
our RG analysis based on the e expansion may shed light 
on this aspect. In particular, unlike some of previous 
works^ii^, the present approach does not rely on the scal- 
ing hypothesis. In stead, the universality in the unitary 
regime naturally arises from the presence of a nontrivial 
zero-density fixed point. Moreover, the roles of danger- 
ously irrelevant operators, like the Ua terms in the action, 
are more easily captured within the RG framework. 



A striking result following from our RG analysis is that 
the T = thermodynamics of this system in the unitary 
limit is identical to that for a free Fermi gas (with a dif- 
ferent mass). This point was also pointed out in Refs. 
and from the variational wavefunction approach. Es- 
pecially, the chemical potential at zero temperature is of 
the form /i = ar?l^ jm in the unitary limit, where a takes 
a universal value. The values of a given by Rcfs. i audi 
are 22.22 and 6.077, respectively. Since the latter is much 
smaller than the former, the authors of Ref. claimed 
that the trial wavefunction adopted by them should be 
an energetically better candidate for the ground state. 
By setting d = 3 in Eq. dSO]), we obtain a = 4.996. 
The fact that the obtained numerical value is close to 
the one in Ref. provides an evidence which further 
supports our approximate solutions based on the e ex- 
pansion. However, to determine the exact value of a, a 
more comprehensive numerical calculation is warranted. 
Moreover, it deserves further study whether or not such 
a correspondence between the dilute FR Bose gas in the 
unitary limit and the free Fermi gas can be extended to 
finite temperature and non-equilibrium processes. 

The quantum phase transition between the ASF and 
MSF phases predicted by the mean-field theorjti can be 
observed only when these phases are thermodynamically 
stable. According to the present analysis (Fig. [1]), the 
ASF and MSF states can be stabilized cither by strong 
dimer-dimer repulsions or by increasing the density of 
atoms to a moderate value. Therefore, to observe such a 
quantum phase transition, we must go beyond the dilute 
limit. 

Our results are obtained by neglecting the higher-order 
correlations such as the three-body effcctsi^. Therefore, 
they are applicable before the Efimov physics fully sets 
^- ^18, 19 ^ It is interesting to see how the three-body effects 
affect the present predictions. 
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